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We point out that, contrary to existing claims of the opposite, Vainshtein's theorem on the non- 
renormalization of the correlator of an axial and two vector currents, is only valid in the chiral 
limit. When quarks, that contribute to the correlator, are massive, the QCD corrections do not 
vanish and the transversal and longitudinal functions are renormalized differently. We compute 
those corrections and study their implications for QCD effects in electroweak corrections to the 
muon anomalous magnetic moment. 



I. INTRODUCTION 

A correlator of an axial and two vector currents {VVA) 
is a peculiar object. Because of its relation to the 
anomaly of the axial current, it plays important role in 
understanding subtle issues in quantum field theory and 
in certain phenomenological applications. 

One such application was recently identified in the 
physics of the muon anomalous magnetic moment where 
calculation of two-loop electroweak corrections requires 
the correlator {VVA), Fig^ Thanks to Furry theorem, 
the vector part of the Z coupling to fermions does not 
contribute and the diagram Figgis entirely determined 
by the {VVA) correlator. Because of significant interest 
in the physics of the muon anomalous magnetic moment, 
driven by recent measurements [l| of this observable at 
Brookhaven National Laboratory, it is not surprising that 
the two-loop electroweak corrections to the muon mag- 
netic anomaly were studied thoroughly. In particular, the 
influence of hadronic effects on electroweak corrections to 
the muon anomalous magnetic moment was discussed in 
a number of recent papers [2, Q| • While these activities 
resulted in a better understanding of hadronic effects in 
the physics of the muon anomalous magnetic moment, 
they also lead to a discovery of an interesting theoreti- 
cal result - a novel non-renormalization theorem for the 
{VVA) correlator 3 (see also @). 

The content of this theorem can be explained as fol- 
lows. In a special kinematic configuration when one of 
the vector currents is soft, the correlator of an axial 
and two vector currents is described by the longitudi- 
nal and transversal functions wl,t- In the chiral limit, 
Wl does not receive QCD corrections as a consequence of 
the Adler-Bardeen theorem on the non-renormalization 
of the axial anomaly [g. The new result, discovered by 
Vainshtein 4], is that in the chiral limit the transversal 
function is also not renormalized by perturbative QCD 



effects. The proof is based on the observation that in the 
chiral limit wl = 2wt, to all orders in the strong coupling 
constant a^. Since the Adler-Bardeen theorem protects 
wl from QCD corrections, the transversal function wt is 
also not renormalized. 

We note that these non-renormalization theorems |j,|a| 
are formulated in the chiral limit. Recently, Pasechnik 
and Teryaev claimed [71] that these results are more gen- 
eral and that even for massive quarks 0(as) QCD correc- 
tions to the longitudinal and transversal functions are ab- 
sent. This result looks somewhat peculiar; for example, 
it is difficult to reconcile the absence of ©(a^) corrections 
with other well-known properties [3, jj| of the correlator 
{VV A) , in particular its operator product expansion. 
Intrigued by that, we decided to re-calculate the QCD 
corrections to the {VVA) correlator in the limit when 
one of the vector currents is soft, allowing for non-zero 
quark masses. In contrast to J3|, we obtain non- vanishing 
QCD corrections away from the chiral limit. Below the 
details of the calculation are described. 



II. CALCULATION 

Consider the correlator of an axial and two vector cur- 
rents 



rpq 

-^ fj.-fi' 



where 



d4xdV^"-^'^(0|T{j^(x)jV(2/)j^(0)}|0), 



(1) 



Jt, = QqQlf^q, jI = 97^759 (2) 

are the vector and axial currents and Qq is the electric 
charge of the quark q. The pole mass of the quark is 
denoted by niq in what follows. 

We assume that momentum k is much smaller than 



q, rUq and view the current j'y{y) as a source of soft pho- 
tons. Then, we can wrfte 



n, = I [d*xe^'^^0\T{j,{x)jtm\^ik)), 



(3) 



where e"'{k) is the polarization vector of a photon with 
momentum k. 




FIG. 1: A two-loop diagram that contributes to electroweak 
corrections to the muon anomalous magnetic moment and 
involves the anomalous {VVA) correlator. 

We are interested in the expansion of the tensor T/^^, 
in k. Such an expansion was constructed in |3,|jl. There 
it was shown that, through first order in fc, T^^ can be 
written 



rpq 



—I 
4^ 



wl{q'^)q„q'' Uf, + w^^iq"^) [-q^f^j,v 



(4) 



qf.q'^ fau - q^q"^ fc 



an 



where /^^ = e^i^apk'^e'^ . The functions w^ rp parametrize 
Lorentz structures contributing to T' , that are longitu- 
dinal (transversal) with respect to momentum q. 






FIG. 2: Examples of two-loop diagrams that lead to 0(0:8) 
correction to w1 ^. The rightmost diagram represents an in- 
sertion of the mass counter-term. 

In general, w^ t depend on q^ and the quark mass 
rUq. In the one-loop approximation, these functions were 
computed in [3,l3 
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a^^OQ^ + ^^l 






(5) 



In Eq.©, iVc = 3 is the number of colors. We also in- 
troduced the Euclidean momentum Q, Q^ = —q^ and 



/3 = J \ + /^w? I Q"^ . In the limit m^ = 0, we find 



w^° = 2w^° = 2iVcQ2/Q2. j^g follows from the Adler- 
Bardeen and Vainshtein's theorems jj, Q , this is exact 
perturbative result in that all QCD corrections to it van- 
ish. 

In this paper, we study 0{as) QCD corrections to Wlx 
in case when the quark mass is not zero. We have to 
compute twelve two-loop diagrams; some examples are 
shown in FigEl Because computing u;' x beyond the 
one-loop approximation for arbitrary value of mi/q^ is 
technically challenging, we adopt a different approach. It 
turns out that the calculation is drastically simplified in 
two limiting cases Q <C rUq and Q 3> fUq. In addition, if 
sufficiently many terms in the expansion of w^ rp in these 
limits are known, we can reconstruct the two functions 
completely using Pade approximation. 



For both cases, Q ^ mq and Q ^ 



^gj 



we employ 



the method of asymptotic expansions ilQj . In the large- 
mass limit, TUq 3> Q, the computation is particularly 
simple since Taylor expansion of Feynman diagrams in 
Fig|2in both k and Q suffices. Consequently, only two- 
loop vacuum Feynman integrals have to be evaluated in 
this case. 

The small mass limit niq ^ Q is more complex, since 
Taylor expansion in mq is insufficient. This can be seen 
from the one loop result which is non-analytic in the limit 
TTiq -^ 0. To expand the two-loop diagrams, Fig[21 in iriq, 
we have to consider how a loop momentum I, flowing 
along a line in a diagram, compares with mq and Q. If 
I ^ Q, the propagator is Taylor expanded in niq. In the 
opposite case, I ~ m^, the propagator l/((g + ^)^ — m^) is 
expanded both in I and rUq while the propagator 1/{P — 
m^) is left unexpanded. 

For each diagram, we have to consider different mo- 
menta routings and identify all possible subgraphs. 
Among them, there are two types of subgraphs that can 
be easily described. The first one corresponds to the sit- 
uation when all lines in a diagram are off-shell by an 
amount of order Q. In this case, Taylor expanding di- 
agram in rriq leads to a massless computation. The op- 
posite situation occurs when both loop momenta are of 
order mq. Then, a diagram is expanded in 1/Q and the 
integrals we have to deal with are the two-loop massive 
vacuum integrals. In intermediate cases when some of 
the lines in a diagram are hard ~ Q and other are soft 
~ TUq, the two- loop graph factorizes into the product of 
one-loop graphs. 

To compute the transversal and longitudinal functions 
separately, we choose different momenta of the soft pho- 
ton. For example, if we choose the soft momentum k in 
such a way that k (x q, we project out the longitudinal 
structure function. To get rid of free Lorentz indices in 
that case, we contract the tensor T^^jy, defined in Eq.©, 



with e'^'^'^'^c 



This procedure allows us to deal 



with scalar Feynman two-loop integrals for the computa- 
tion of w^. As we pointed out in the previous paragraph, 
the calculation of these integrals is simple once the limits 



rUq <^ Q OT niq ^ Q are considered. A similar procedure 
is used to compute w^. 

We use dimensional regularization for the computa- 
tion. For the consistent treatment of the axial current in 
d dimensions, we employ the procedure of Ref. llj. The 
calculation is performed in an arbitrary covariant gauge; 
the cancellation of the gauge parameter dependence in 
the final result is a welcome check on the correctness of 
the calculation. To present the results for the longitudi- 
nal and transversal functions, we define 
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For Q <^ niq, we find 



A 



9:(1) 



A, 



9,(1) 






617 Q^ 8041 Q^ 



540 m'^ 



25200 m4 



521 Q2 449 Qi 



540 to2 



1680 m4 



(7) 



where ellipses stands for higher order terms in the ex- 
pansion in Q^ /mi. In case Q 3> iriq, the functions read 
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where L — ln{Q'^/m'i), C(3) is the Riemann zeta-function 
and ellipses stands for other terms suppressed by higher 
power of m'^/Q^. 

Several things can be pointed out in connection with 
these results. As expected, the QCD corrections to w^ j, 
vanish in the limit rriq — > 0, in accord with the non- 
renormalization theorems |j,|g. Also, the QCD correc- 
tions are not universal so that the all-orders perturba- 
tive relation w'[ — 2w^ is violated, once non-zero quark 
masses are allowed. 

Using Eas. (|7l8|l . we can construct the entire functions 
A^ T with the help of Fade approximation 12j. To this 
end, we compute ten terms in the expansion of A^ t in 
Q^/m^ and perform Fade approximation replacing the 
series by rational functions. It turns out that the re- 
sulting functions do not depend on the precise form of 
the Fade approximation employed for their construction. 
For computations described below we use the [4/5] Fade 
approximant. 

The virtue of the Fade approximation is that it allows 
us to continue the small-Q expansion of the functions 



FIG. 3: The functions A^ (upper curve) and A^ (lower 
curve) plotted in dependence on a; = Arriq/Q^. 



A^^x to large values of Q, where it should smoothly 
merge with the large-Q expansions of these functions 
Eq.®. Indeed, we observe a perfect match of the 
Fade approximants and the large-Q^ expansion of the 
functions A^ t i^^ ^^^ interval lOm^ < Q^ < lOOm^. 
Hence, the Fade approximants, supplemented with large- 
Q^ asymptotic, can be used to deduce A^ t ^^^ arbitrary 
values of Q^/m^. 

For completeness, we present below the functions 



,9,(1) 



in the form of the [4/5] Fade approximants. We 



introduce z = Q /m and write 



A 



9,(1) _ 
L,T — 



4 
i=l 



(*) , 



J2 ^L,T^ 
i=l 



(9) 



The coefficients aj^'^, , bi^'j, are given in Tables HITll The 
functions A^ x a-ie plotted in Fig|3| It is interesting 
that for large range of m^/Q^, the approximate relation 



,9,(1) 



2.32AII;' ' holds with reasonable accuracy. 





L 


T 


a(i) 


1.1666666666 


0.5000000000 


a(^' 


0.3762074740 


0.1620591351 


a(3) 


0.03273978142 


0.01418000768 


a(4) 


6.50545897 ■ lO"* 


2.836663049 • 10~* 



TABLE I: Coefficients for the [4/5] Fade approximant for the 
longitudinal and transversal functions. 



III. RELATION TO THE OPE 



As we pointed out in the Introduction, the non- 



renormalization of the functions 



^L,T 



by perturbative 
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T 


&(1) 


0.6489185756 


0.6457232084 


6(2) 


0.1487376001 


0.1469404864 


5(3) 


0.01420626517 


0.01388076895 


5(4) 


5.028490259 ■ 10"" 


4.828123552 • 10~* 


5(5) 


3.78544 • 10"'' 


3.527956921 ■ lO"*^ 



TABLE II: CoefRcients for the [4/5] Pade approximant for 
the longitudinal and transversal functions. 



to compute quadratic mass correction to the massless 
limit of the transversal and longitudinal functions. The 
first one is the dimension-two operator, O^^ = Fa/3/(47r), 

where Fap = e-aPfivd^A'^] its matrix element between 
the vacuum and the soft photon is trivial and leads to 
Kp ~ 1. The second operator is the dimension-three 
operator O'^ = —iqcTafSjzq- Its leading order Wilson 
coefficient reads 0, Q 



cl{f,) = 24ip) = 



4:Qqmqip) 



(12) 



QCD effects beyond the chiral limit, claimed in Q, is 
not possible to understand given the operator product 
expansion (OPE) of the tensor T/^^ derived in 0,3- In 
this Section, we elaborate on this statement. 

Typically, the OPE is employed to estimate non- 
perturbative corrections to, otherwise, nearly pcrturba- 
tive observables. This is achieved by separating physics 
at different distance or momentum scales. Soft contribu- 
tions, associated with non-perturbative effects, are char- 
acterized by momentum scale Aqcd whereas hard con- 
tributions are characterized by external kinematic scale 
Q » Aqcd- When the OPE is applied, soft compo- 
nents arc identified with matrix elements of local oper- 
ators while hard momenta contribute to Wilson coeffi- 
cients of these operators. 

Although the OPE is associated with computations be- 
yond perturbation theory, its main idea, separation of 
physics at different distance scales, makes it suitable for 
perturbative computations if largely different scales are 
present. Hence, it is well-suited for studying longitudi- 
nal and transversal functions in the limit Q ^ iriq. In 
this case, soft scale is set by the mass of the quark m^ 
rather than Aqcd, but the OPE remains intact. 

The OPE of the longitudinal and transversal functions 
reads @,[| 



where mq(fi) is the running quark mass 



f'L.T^ 



= y cl^T.{p)K,{p), 



(10) 



A^^)^ 



1 + 



Cpas 



3 M^ 

-m — TT 



(13) 



From the OPE of the functions wl ^, Eq.l^UI), it fol- 
lows that, in order to compute these functions through 
0{asm'i/Q'^), we require the Wilson coefficient of the 
operator Ofa, the matrix element of the operator O^^ 
and the Wilson coefficients c^ -j through 0{as)- To ob- 
tain those ingredients, two- and one- loop computations 
are needed. It turns out that all the matrix elements 
and Wilson coefficients can be extracted from the results 
reported in the previous Section. This happens because 
our calculation is based on computing contributions from 
different momenta regions separately. Then, soft regions 
in the one- and two-loop cases are identified with the 
matrix element of the operator O'^n through 0{as). We 
have checked that identical results for this matrix element 
are obtained, independent of whether such identification 
is done in wf^ or w!^. This independence is simple, yet 
welcome, check on the correctness of the calculation. 



/ 




where the normalization scale fi is introduced to make 
the separation of hard and soft modes unambiguous. In 
Ea. (|10|l . the sum includes all local operators that con- 
tribute to the matrix element (0|C'*|7(fc)) and behave as 
rank two pseudo-tensors under Lorentz transformations. 
The matrix elements are parametrized by 



\oi,f,{pMk)) 



-^faf3lii(p), 



(11) 



FIG. 4: The one-loop contribution to the matrix element 
of (0|C'^g|7(fc)). The two-loop contribution is obtained by 
adding all possible single gluon exchanges to this diagram. 



The matrix element (010^017) is divergent (cf. Fig0J. 
To arrive at the renormalized matrix element we use 
the MS renormalization scheme and the fact that, under 
renormalization, O' „ mixes with the operator mFap. We 
obtain 



where the dependence on the normalization point is made 
explicit. 

From our results in the previous Section, it follows that 
the relation wl = 2i(;t is violated once 0{asiTi^ /Q^) con- 
tributions to these functions are computed. We would 
like to analyze this result from the OPE perspective. Ac- 
cording to |3,y| , there are two operators that are required 



Kq{fl) = NcQqmq{fJ,) 


8 


Cpasifi) ( Ll + 

IT \ 4 


Lf, 119 TT^ 

+ 48 + 8 



where L^ = ln{fi^ /ni^di)). 



(14) 



Using Eq. (|14|l , we determine Wilson coefficients of the 
operators O^'S by matching the OPE of the functions 
, Eg. lfTUIl . to perturbative result Eq.®. For simplic- 



%,T 



ity, we give those Wilson coefficients at the normalization 
scale n = Q, where they do not contain any logarithms. 
To present the results, we introduce a short-hand nota- 
tion fhq — mq{Q) and 



Cl,t 






(15) 



Through 0{as), wc find 



c[-2-^' 



2 

2 r 



c|: = i 



and 



Q' 



32 Cpcts 
16 Cpcts 



4:Qqmq 
2Qqfhq 



-SCa 



61 

24 



T 



4Cc 



(16) 



Q' 



CpOs 

27r 



(17) 



We use TUc = 1.5 GeV, nib = 4.8 GeV and mt = 
180 GeV for charm, bottom and top quark masses. We 
also employ the leading order running of the strong cou- 
pling constant with as^rric) = 0.3. We obtain the fol- 
lowing QCD corrections due to charm, bottom and top 
quarks 



a Gum,, 

'' ttStt^V2 



1.80, 

-0.33, 

0.56, 



q = c, 
q = b, 
q = t. 



(20) 



The perturbative QCD correction is given by the sum of 
entries in Ea. lf^ . We obtain 



Aa"= = 0.55 X 10" 



(21) 



This correction is well within the error bars, assigned 
to hadronic uncertainties in electroweak contributions to 
the muon magnetic anomaly in Ref . J3] ; it is negligible for 
phenomenology given that the current uncertainty on the 
muon anomalous magnetic moment is ~ 100 x 10~^^. 



CONCLUSION 



IV. THE MUON ANOMALOUS MAGNETIC 
MOMENT 

The correlator of an axial and two vector currents in 
the kinematic limit when one of the vector currents is 
soft, appears in the physics of the muon anomalous mag- 
netic moment, Fig^ The correction to the muon mag- 
netic anomaly due to these diagrams is obtained 2^ upon 
integrating w'[ rp over Q^ with the weight function 



Aa„ 



_ aGjjn^ 
''^~ 7r87r2y2 

In Ea. (fTH|) we use 



dQ" wi^ 



^WT 



(18) 



WL.T 



-E2^: 



/ / 

3'"'l,T' 



(19) 



and the summation index / denotes both lepton and 
quark contributions to the (VVA) correlator. In addi- 
tion, /g stands for the weak isospin of the fermion /. 

In Section^ we derived Wlt' including 0[as) QCD 
corrections. We can use those results in Ea. lfTH|l to ob- 
tain QCD corrections to Aa^. Note that, at large Q, 

the functions A^ t decrease as m^/Q^; hence, at order 
0{as) we can compute the integral in Ea. H18l) for each 
quark fiavor separately. We only include charm, bottom 
and top quarks in the computation. The up, down and 
strange quarks have to be treated non-perturbatively [2| ■ 



In this paper we computed the QCD corrections to 
the longitudinal and transversal functions w'^ rj, allowing 
for non-zero quark masses. These functions describe the 
anomalous {VVA) correlator in the kinematic limit when 
one of the vector currents is soft. 

According to the Adler-Bardeen and Vainshtein's the- 
orems, Wl t do not receive QCD corrections in the chiral 
ruq = limit. Our study is motivated by the claim in 
Ref. PI that the absence of QCD corrections w^ x is also 
valid beyond the chiral limit when quark masses are al- 
lowed. 

We performed explicit calculation of the QCD cor- 
rections to 'w'[rj,{Q) for non-zero quark masses in two 
kinematic limits niq <C Q and Q <C mq and used Fade 
approximation to derive the approximate form of these 
functions valid for arbitrary relation between Q and niq. 
In variance with results of Ref. |3] , we found that w'[ rp 
receive perturbative corrections if niq ^ 0. We ar- 
gued that these QCD corrections are consistent with the 
well-known operator product expansion of the correlator 
iyV A). We studied the implications of our results for 
the QCD corrections to electroweak contributions to the 
muon anomalous magnetic moment that contain the cor- 
relator of an axial current and two vector currents, Fig^ 
and found negligible effect. 
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